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I. GRIBOV LIPATOV RECIPROCITY 

Evolution equations describing dynamics of 
distributions of QCD partons (quarks and gluons) 
that determine cross sections of hard processes 
have a long history. 

Logarithmic violation of the Bjorken scaling 
- dependence of DIS structure functions on the 
hard scale of the process - and, simultaneously, 
scale dependence of fragmentation functions in 
e + e~ annihilation, were first addressed in a gen- 
eral QFT context by Gribov and Lipatov back in 
19710. 

These fundamental results were explicitly 
presented in the language of evolution of a system 
of partons by Lipatov in 1974 [2]. QCD parton 
evolution equations followed suit in 1977 [3,13]. 

The universal nature of the parton dynamics 
goes under the name of factorisation of collinear 
("mass") singularities. 

Physically, it is due to the fact that quark- 
gluon multiplication processes happen at much 
larger space-time distances than the hard inter- 
action itself. It is this separation that makes 
it possible to describe quark-gluon cascades in 
terms of independent parton splitting processes. 
They success one another in a cleverly chosen 
evolution time, t ~ InQ 2 , whose flow "counts" 
basic parton splittings that occur at well sepa- 
rated, strongly ordered, space-time scales. 

Perturbative structure of the cross section of 
a given process p characterised by the hardness 
scale Q 2 can be cast, symbolically, as a prod- 
uct (convolution) of three factors (for a review 
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see @|): 



= cW[a,(t)]®exp(jf drP[a s (T)] 
<8> Wh(to), t ~ InQ 2 . 
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Here the functions C[a s ] (hard cross section; co- 
efficient function) and P[a s ] (parton evolution; 
anomalous dimension matrix) are perturbative 
objects analysed in terms of the a s -expansion. 
The last factor Wh embeds non-perturbative in- 
formation about parton structure of the partic- 
ipating hadron(s) h, be it a target hadron in 
the initial state (parton distribution) or a hadron 
triggered in the final state (fragmentation func- 
tion). 

Gribov and Lipatov found the splitting func- 
tions Pab (x) determining the evolution of parton 
fluctuations to be identical for space-like (DIS) 
and time- like (e + e~ annihilation) cases: 
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The identity of "evolution Hamiltonians" in the 
two channels is a highly non-trivial property since 
the "energy fraction" arguments the splitting 
functions depend on ([3]) are not the same but 
rather reciprocal. 

Combined together with the Drell-Levy- 
Yan relation [|| (which followed from similarity 
of Feynman diagrams for the two crossing chan- 
nels), the identity ([2]) translated into an internal 
symmetry relation for the splitting functions, 

(T)x ■ Pba{x~ 1 ), 



Pab(x) 
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known as the Gribov-Lipatov reciprocity (with 
the =f 1 factor depending on the spins of partici- 
pating partons A, B; for a review see, e.g. 0]). 
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However, the Gribov-Lipatov reciprocity 
(GLR) was found to be broken beyond the leading 
logarithmic approximation (the one-loop parton 
Hamiltonian) . 

Having observed the difference between 
spase- and time-like two- loop splitting functions, 
Curci, Furmanski and Petronzio [8j have re- 
marked that the violation of the GLR might be 
of kinematical origin. 



II. RECIPROCITY RESPECTING 
PARTON EVOLUTION EQUATION 

A few years ago the follo wing evolution 
equation has been suggested in 0, 1 1 Of ] : 



dD(x,Q 2 ) 
dlnQ 2 



J ^V[z,a s ]-D(- 2 ,z°Q 2 
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where a = — 1(+1) for the space-like (time-like) 
case. Appearance of the combination z a Q 2 in 
the hardness argument of the parton distribution 
under integral corresponds to choosing fluctua- 
tion lifetimes of successive virtual parton states 
as a logarithmic ordering parameter ("evolution 
time"). 

Such a choice of an evolution variable is 
known to be "wrong" for either of the two chan- 
nels. 

Instead, it is the transverse momentum or- 
dering, k 2 ^ in the S channel, and the angular or- 
dering, (k±/z) 2 , for the T-evolution, correspond- 
ingly, that make anomalous dimensions free of 
series of double logarithmically enhanced terms 
(a s \nx 2 ) n . Such terms become explosively large 
in the small- x region which region is practically 
important both for DIS (high energy scattering) 
and for jet physics (soft gluon multiplication). 

At the same time, the fluctuation time or- 
dering variable, k\jz, happens to be lying right 
in between the two "clever choices" and thus pre- 
serves the symmetry between the two channels. 
Therefore, one may hope that formulating par- 
ton evolution in terms of equation ([5]) will result 
in the universal (channel independent) reciprocity 
respecting evolution kernel V(x). 

Non-locality of the new equation ([SD in lon- 
gitudinal (z) and transverse variables (hardness 
scale Q 2 ) breaks identification of splitting func- 
tions with anomalous dimensions. 



What it offers instead is a link between the 
two channels by means of universal evolution ker- 
nel matrix V, one and the same for T and S evo- 
lution. In spite of the fact that the new "split- 
ting functions" V in (|5|) do not correspond to any 
clever choice of the evolution variable, in either 
T or S channel (explosive a s In 2 x terms being 
present in both cases), this universality can be 
exploited for relating DIS and e + e~ anomalous 
dimensions. 

One can expect that by separating the no- 
tions of splitting functions and anomalous di- 
mensions by means of the Reciprocity Respect- 
ing Evolution Equation ([5]) the Gribov-Lipatov 
wisdom can be rescued in all orders. 

This guess was motivated by a remark 
made by Curci, Furmanski & Petronzio @| who 
observed that the GLR violation in the sec- 
ond loop non-singlet quark anomalous dimension 
amounted to a "quasi- Abelian" term oc C% with 
a suggestive structure 



P${z)hiz. 
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Their observation hinted that the GLR violation 
was not a dynamical higher order effect but was 
inherited from the previous loop via a non-linear 
relation [9|]. 

In the Mellin space the convolution ([6|) 
translates into 



PnP, 



N *N • 



Let us check that it is this structure of the GLR 
breaking that emerges from @. 

Taking Mellin moments of both sides of the 
equation we obtain 



la (Np N (Q 2 



1 A~ 

-z N P[z,a s ]z^D N (Q 2 ), 



where we have used the Taylor expansion trick. 
The integral formally equals 

la {N) = {D^ViN + ad^DN, (7) 

expressing the anomalous dimension through the 
Mellin image of the evolution kernel with the dif- 
ferential operator for the argument: 

N^N + ad lnQ2 . 
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The derivative acts upon Dn{Q 2 ) producing, by 
definition, j(N)D^. In high orders it will also 
act on the running coupling the anomalous di- 
mension depends on, 7 = i(N, a s ). The latter 
action gives rise to terms proportional to the (3- 
function. Such terms are scheme dependent as 
they can be reshuffled between the anomalous di- 
mension and the coefficient function C[a s ] in the 
expression ([T]). 

Neglecting for the time being such contribu- 
tions by treating a s as constant, ([7]) reduces to a 
functional equation 



7ct (A0 = V{N + a la {N)) 



(8) 



Since 7 = 0(a s ), we can expand the argument of 
the evolution kernel perturbatively, 

7<r = V+V-a 1 +\V- 1 2 + 0{l3{a)) + 0{a i ). (9a) 

Solving ([9a]) iteratively we get 



la 



V + aVV + 



VV 2 + \V 2 V 



+ ... (9b) 



Restricting ourselves to the first loop, V = a , 
with the (Mellin image of) good old LLA 
functions, gives 



la 



a 



pW + aPapWpV) + ... (10) 



The second term on the r.h.s. of (|10p generates 
the two-loop Curci-Furmanski-Petronzio relation 
© all right. 

Knowing the n-loop anomalous dimension 
matrix in the S channel, the RREE predicts the 
anomalous dimensions in the T channel (and vice 
versa) . 

Based on the existing three-loop results for 
the space-like evolution [111 ], the corresponding 



prediction of © for the time-like channel was 
verified in the case of non-singlet anomalous di- 
mensions by Mitov, Moch and Vermaseren in [T3 |. 

Basso and Korchemsky have traced the ori- 
gin of the relation ([8]) to the underlying confor- 
mal properties of the theory and described how 
to embed into the equation effects of non-zero (3- 
function [T3 |. 

Validity of the Gribov-Lipatov reciprocity 
for the kernel V of the new evolution equation 
([5]) has been checked not only for three-loop non- 
singlet QCD anomalous dimensions. 

Basso and Korchemsky (in collaboration 
with Moch) have revisited virtually all known 
multi-loop QFT results and found GLR to hold 
for three-loop unpolarized singlet and two-loop 
polarized QCD distributions (quark transversity, 
linearly polarized gluon, quark singlet polarized), 
X(j) QFT at four loops, QCD in the /3o — > 00 
approximation (in all loops), as well in vari- 
ous SUSY models, including multi-loop anoma- 
lous dimensions of the maximally supersymmetric 
M = A YM theory. In the latter model the GLR 
was found to hold even in the strong coupling 
limit, a — > 00 (accessible through the AdS/CFT 
correspondence) . 



III. QUASI-ELASTIC LIMIT 

The RREE's first demonstration of force 
was the derivation of all-order predictions for 
the structure of subleading singular terms in the 
expansion of the quark non-singlet and gluon 
anomalous dimension in the large- x limit [101 ] . 

Behaviour of anomalous dimensions in the 
quasi-elastic (x — > 1) kinematics can be cast in 
the following form: 



la(x) 



Ax 



(1-x) 



+ B5(l-x) +C a ln(l-x) + D a + 0({l-x)\og p (l-x)) 
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where the coefficients A and B are the same in singular — term x/(l — x) is dictated by the cel- 
the two channels. ebrated Low-Burnett-Kroll (LBK) theorem 14]. 

Specific structure of the first — the most It is a consequence of the fact that soft radiation 
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at the level of da oc duoiuj^ + const) has classical 
nature. 

The coefficient A in front of this structure 
has a meaning of the "physical coupling" mea- 
sured by the intensity of relatively soft gluon 
emission. This coefficient (calculated in three 
loops in the MS-bar scheme) is known to uni- 
versally appear in all observables sensitive to soft 
gluon radiation: quark and gluon Sudakov form 
factors and Regge trajectories, threshold resum- 
mations, singular part of the Drell-Yan iT-factor, 
distributions of jet event shapes in the near-to- 
two-jet kinematics, heavy quark fragmentation 
functions, etc. The structure (fTTj) applies to the 
large- x behaviour of the g — > g anomalous di- 
mension as well, with Ar g \/At q \ = Ca/Cf, in all 
orders. 

Quantum effects show up only at the level 
of da oc ijj duj that is at the level of contribu- 
tions that were neglected in (fTTj) . At one loop, 
subleading terms Cln(l — x) and the constant D 
in (jlip are absent. This suggests that in higher 
loops they should emerge as "inherited" rather 
than non-trivial entries. 

Indeed, to keep under control all the terms 
in ([i~l~]) it suffices to use the x — > 1 asymptote 
of the one-loop evolution kernel = to 
derive from RREE ([5]) all-order relations [10f| 



-aA 2 , 
-a AB. 



(12a) 
(12b) 



The relation (fT2"a|) is "conformal" while (pbj) 
acquires correction due to running of the cou- 
pling [la ]. 

The ideas of the universal evolution equation 
have recently found an interesting application in 
the work by Laenen, Magnea and Stavenga who 
have employed the RREE as means of improving 
threshold resummations [l| 



IV. Af = 4 SYM 

QCD shares the gluon sector with supersym- 
metric YangMills models (SYM). This suggests to 
explore supersymmetric partners of QCD in order 
to shed light on the subtle structure of the per- 
turbative quarkgluon dynamics. 

QCD is not an integrable quantum field the- 
ory. In spite of this, in certain sectors of the chro- 



modynamics the integrability does emerge [R 
This happens, markedly, in the problem of high 
energy Regge behaviour of scattering amplitudes 
in the large-iV c approximation (planar 't Hooft 
limit), in the spin 3/2 baryon wave function, for 
the scale dependence of specific (maximal helic- 
ity) quasi-partonic operators (for review see [TtJ )• 
What all these problems have in common, is the 
irrelevance of quark degrees of freedom and the 
dominance of the classical part of gluon dynam- 
ics, in the sense of the LBK theorem [lit]. 

The higher the symmetry, the deeper inte- 
grability. The maximally supersymmetric Af = 4 
YM theory is exceptional in this respect. A st ring 
of recent theoretical developments 13, E, Gif 
hinted at an intriguing possibility that this QFT, 
super-conformally invariant at the quantum level 
(/3(a) = 0), may admit an all loop solution for 
anomalous dimensions of its composite operators. 

The Af = 4 SYM being an integrable model, 
there exists a powerful technology based on the 
Bethe Ansatz Equations well suited for perturba- 
tive calculation of anomalous dimensions of com- 
posite operators to multi-loop accuracy. The so 
called "universal anomalous dimension" of the 
Af = 4 SYM theory is given by the "maximal 
transcedentality" Euler-Zagier harmonic sums 

Applied to this theory, the kernel of ([5]) was 
found to respect GLR in four loops for the leading 
twist two [221 , as well as in four (23|] and five loops 



241 ] for twist three operators. 

The RREE was found to significantly sim- 
plify the structure of high order terms in the "uni- 



versal anomalous dimension" of the theory 251 ] . 

Given that in the leading order evolution 
kernel of the Af = 4 SYM is purely classical in 
the LBK sense, 



x 



- + no quantum corrections, (13) 



one may hope to derive one day a one-line-all- 
loops expression for the anomalous dimension of 
this theory, in which higher order terms are dy- 
namically "inherited" from the first loop. 

QCD would greatly benefit from such a so- 
lution, since this Af = 4 SYM finding would put 
under full theoretical control the dominant part 
of the perturbative QCD gluon dynamics. 
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